Observation of non-exponential magnetic penetration profiles in the Meissner state — 
A manifestation of non-local effects in superconductors 
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Implanting fully polarized low energy muons on the nanometer scale beneath the surface of a 
superconductor in the Meissner state enabled us to probe the evanescent magnetic field profile B(z) 
(0 < z < 200nm measured from the surface). All the investigated samples [Nb: k ~ 0.7(2), Pb:« ~ 
0.6(1), Ta: n ~ 0.5(2)] show clear deviations from the simple exponential B(z) expected in the 
London limit, thus revealing the non-local response of these superconductors. From a quantitative 
analysis within the Pippard and BCS models the London penetration depth Al is extracted. In 
the case of Pb also the clean limit coherence length £o is obtained. Furthermore we find that the 
temperature dependence of the magnetic penetration depth follows closely the two-fluid expectation 
1/A 2 oc I — (T/T c ) 4 . While B(z) for Nb and Pb are rather well described within the Pippard 
and BCS models, for Ta this is only true to a lesser degree. We attribute this discrepancy to the 
fact that the superfiuid density is decreased by approaching the surface on a length scale £o- This 
effect, which is not taken self-consistently into account in the mentioned models, should be more 
pronounced in the lowest k regime consistently with our findings. 

PACS numbers: 76.75.-K, 74.20.-z, 74.25. Ha, 74.25.Nf, 74.78.-w 
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I. INTRODUCTION 



One of the fundamental properties of a superconductor 
is the Meissncr-Ochsenfeld effecti It states that at low 
magnetic fields and frequencies a superconductor expels 
or excludes any magnetic flux from its core. However, at 
the surface the field penetrates on a typical length scale 
A called the magnetic penetration depth. In the Lon- 
don limits, for a semi-infinite superconductor, the func- 
tional B vs. z dependence is an exponential one. This 
approximation holds for a large class of superconductors. 
However, as found by Pippard from a set of microwave 
experiments on impurity doped type-I superconductors^, 
this is not always true. In analogy with the anomalous 
skin effect Pippard introduced the concept of non-local 
response of the superconductor, i.e. the screening cur- 
rent trying to expel the magnetic field must be averaged 
over some spatial region of the order of £ called the co- 
herence length. The physical interpretation of £ is, that 
it is the length over which the superconducting wave 
function can be considered as rigid, i.e. roughly speak- 
ing the size of a Cooper pair (for details see Sec. ITT|) . 
The non-local electrodynamical response leads to vari- 
ous modifications of the London theory, one being that 
the magnetic penetration profile B[z) is no longer expo- 
nential and even changes its sign beneath the surface of 
the superconductor. All these findings were confirmed by 
the microscopic BCS theorjii. Although these theoretical 



predictions have been known for half a century, only very 
recently has a "direct" measurement of the functional 
dependence of B(z) been demonstrated^. A historical 
summary of the experimental work on this subject be- 
gins with that of Sommerhalder and co-workeii&iS who 
showed the existence of a sign reversal of B(z) by mea- 
suring the magnetic field leaking through a very thin su- 
perconducting tin film. Doczcma et alM> applied magnc- 
toabsorption resonance spectroscopy techniques to tackle 
the problem. The technique uses the fact that quasi- 
particles traveling parallel to the shielding current are 
bound to the surface by an effective magnetic potential. 
Indication of non-local effects in Al were inferred by com- 
paring microwave induced resonant transitions between 
the energy levels of these bound states with transition 
fields calculated from the energy levels of the trapping 
potential, parameterized to include the shape of the non- 
local BCS-like potential. Due to the resonant character 
of the experiment, only a few specific points of the po- 
tential are probed. In addition the normal metallic state 
has to be understood very well in order to interpret the 
data. This, together with uncertainties in modeling the 
surface bound states, leaves room for speculations. Po- 
larized neutron reflcctometry has also been applied since 
specular reflectivity of neutrons spin polarized parallel or 
anti-parallel to B depends on the field profile. However, 
this technique requires model-fitting of spin-dependent 
scattering intensities rather than giving a direct measure 
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of the spatial variation of the magnetic field. Up to now 
non-local corrections have been found to lie beyond the 
sensitivity of polarized neutron reflectivity techniques^. 
In this paper we present magnetic field profiles measured 
beneath the surface of various superconductors in the 
Meissner state by means of low energy muon spin rota- 
tion spectroscopy (LE-//SR). The results provide a direct 
and quantifiable measure of non-local effects in the inves- 
tigated materials (Nb, Pb, Ta) and permit the extraction 
of physical parameters such as the magnetic penetration 
depth A and the coherence length £. The paper is orga- 
nized as follows: The next section reviews the theoretical 
framework necessary to understand the results and the 
discussion. Section ITTT1 provides some information on the 
experimental technique including the characterization of 
the samples. In Sec. HVI wc present our data, including a 
discussion, followed by a summary in Sec. [3 

II. THEORETICAL BACKGROUND 

In this section, the theory of the response of a su- 
perconductor in respect to an external electromagnetic 
field will be sketched (for a rigorous derivation see Ref. 
ITlf) . An external electromagnetic field acts on the ground 
state of a superconductor as a perturbation^. Within 
standard perturbation expansion one can showii that the 
following non-local relation between the current density 
j and the vector potential V A A = B, where B is the 
magnetic induction, holds: 

Ur) = E / { R M ~ ^T- <^)<W }Mr') dr' 
P ™ ' 

=: K aP (p) 

(1) 

where p = r — r', e is the charge, ng the supercurrent 
density, and m* the effective mass of the charge carrier. 
K a p{p) is called the integral kernel. The vector potential 
A(r) needs to be properly gauged in order that Eq. Q 
is physically meaningful^. R a p(p) describes the param- 
agnetic response, whereas the second term in the bracket 
reflects the diamagnetic one. If the ground state wave 
function of the superconductor were "rigid" with respect 
to all perturbations (rather than only those which lead 
to transverse excitations) R a p would be identically zero 
and Eq. would reduce to the second London equation 

3a(r) = -—^A a (r) (2) 

Mo A L 

with the London penetration depth Al = 
y m* I (noe 2 ns), and /jq the permeability of the vacuum. 
This, together with the Maxwell equation VAB = fioj, 
results, for a semi-infinite sample, in the well known 



penetration profile 



B(z) = B oxt exp(-z/A L ) (3) 
where z is the depth perpendicular to the surface and 
-B ex t the externally applied magnetic field strength. 

In situations where the paramagnetic term R a p(p) in 
Eq. £Q| cannot be neglected one arrives at the more com- 
plicated formula 

B(z) = B cxt - [ g sm(qz)dq. (4) 

K(q£, T, £) is the Fourier transformed kernel from Eq. JJJ 
including the electron mean free path I. Since only the 
one dimensional case will be considered everything is ex- 
pressed in scalar form. This equation reduces obviously 
to an exponential decay if K(q£,T,£) is independent of 
q, and in the London limit fiQK(q^T,£ — > oo) = 1/A|\ 
Eq.Q is valid in the case of specular reflection of the 
charge carriers at the surface. Another extreme limit is 
given by assuming diffuse scattering at the interface. A 
real system will not be exactly in one of the two limits. 
Since a quantitative analysis^ shows that the difference 
in A is marginal, we use Eq.(@J for the following discus- 
sion. Still, we have verified by numerical integration of 
Eq.Q that B(z) for specular and diffuse scattering are 
closely related. 

Since the electromagnetic response of the supercon- 
ductor is on the length scale of A, a local description 
would be sufficient if the kernel noK(q£, T, £) is constant 
in the interval < q < 1/A. In the next paragraphs we 
will show that ^K{qS > ,T 1 (.) is approximately constant 
for < q < l/£. From this it follows, that for A » £ 
(type II) the local London limit should be quite reason- 
able, whereas in the opposite limit A £ (type I) this 
is not the case. For the type II case, non-local effects 
might play a role cither in case the energy gap exhibits 
nodes and henceforth £ will be very anisotropic^ or if 
K = A/£ < 1.4. We would like to stress that the border 
line between the non-local and the local regime (k w 1.4) 
does not coincide with the border line between type-I and 
type-II superconductivity (re = 1/^/2 = 0.7071 . . .). 

Knowing K(q£,T,£) enables one to calculate B(z). 
The functional dependence of K(q£,T,£) was derived 
scmi-phcnomcnologically by Pippard and later micro- 
scopically by Bardecn, Cooper and Schricffer. 



A. Pippard Kernel K P (q£,T,£) 

Starting from a possible analogy between the anoma- 
lous skin effect and the Meissner-Ochscnfeld effect, 
Pippard^ arrived at the formula for the kernel 
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fioK P (q£,T,e) 



A 2 (T) &(0,£) 



3 1 

2^3 



{ (l + x 2 ) arctan(x) — a;} 



(5) 



with x = q£p (T, t) , and the Pippard coherence length 
£p (T, t) . The temperature dependence of £p (T, could 
be explained only later by the BCS theory (see also the 
next section) and is 
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J(0,T) , 1 



fr(T,t) fr(0) 



(6) 



with 



A(T) 
2k B T 



where A(T) is the superconducting energy gapi* and fee 
the Boltzmann constant. The weak temperature depen- 
dence of £p(T, £ — > oo ) is shown in Fig^ 




Figure 1: Temperature dependence of the Pippard coherence 
length £p(T,£ — > oo) according to Eq. (J^J. The inset shows 
the temperature dependence of J(0, T) within the weak cou- 
pling limit. 



B. BCS Kernel K B cs(q£,,T,£) 



with the following set of abbreviations 



K{T,£) 
1 



In 
A(0) 



UT,£) 



t 



2 fn A(T) 1 
tt Co A(0) + * 



(8) 



A(T) : 



ttA(O) 



v/l + (2n + l) 2 a 2 



k B T c = 1.764 fc B T c 



as 



The temperature dependence of A is defined 
1/A 2 (T) := \im q ^ ^K(qtT,£), i.e. A(T -» 0,£ - 
oo) = Al, and £o the clean limit coherence length at 
T = 0. Though the BCS expression is much more in- 
volved, the q dependence of the kernel is very close to 
the one given in the phenomenological Pippard expres- 
sion Eq. (JSJ. A comparison is given in FigEI 
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Figure 2: (Color online) (/-dependence of the different ker- 
nels for T — > 0. The dashed line represents the g-independent 
London approximation. The BCS- (full) and Pippard-kernel 
(dash-dotted) are very similar and have a q~ asymptotic be- 
havior for q — > oo. 



Starting from the weak coupling BCS model*, one ar- 
rives at the following expression for the kernel 15 



/'n 



OO 

K BCS (qt, T,£) = Y / jjjTfi ■ 9[qUT, t)} (7) 



n=Q 



The corresponding magnetic fields can only be calculated 
numerically. An example of an extreme type-1 super- 
conductor, Al, is presented in Fig|3J As intuitively ex- 
pected compared to the local case, the initial slope of 
B(z) is reduced, reflecting the fact that the magnetic 
field penetrates deeper into the superconductor. The 
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Figure 3: (Color online) The magnetic penetration profile in 
the Meissner state of Aluminum according to Eq.Q for a 
£o = 1600 nm and Al = 50 nm. The inset shows the same 
curve in a log scale. 

log \B(z)/ B cxt \ inset of FigEl shows initially a clear neg- 
ative curvature, i.e. a deviation from the exponential 
behavior. The next important feature to be noticed is 
the sign reversal of the magnetic field before approach- 
ing zero deep inside the sample. All these findings can 
be made plausible by the following hand waving argu- 
ments: In the non-local case the Cooper pairs are very 
extended compared to the magnetic penetration profile. 
Since the partners within a Cooper pair do not experience 
the same field, the screening response is less effective and 
hence the slope is less steep compared to a local response. 
This has a second effect: since the field penetrates further 
beneath the surface, at some depth enough Cooper pairs 
will experience it and start to "overcompensate" , which 
accounts for the negative curvature as well as for the field 
reversal of B(z) before approaching zero. This pictorial 
view is clearly an oversimplification since the screening is 
due to the Cooper pairs themselves, a strongly feedback 
coupled system, but it provides the basis for a qualitative 
understanding of the penetration profile. 

C. Strong Coupling Corrections 

The BCS weak coupling theory assumes a weak 
electron-phonon interaction. This is definitely not the 
case for either Pb or Nb. Therefore, the much more 
elaborated strong-coupling thcorjii& has to be used to 
calculate the properties of the superconductor. A recent 
review describing the strong coupling theory is found in 
Ref. 0. Fortunately, the outcome for the magnetic pen- 
etration profile B(z) in the Meissner state is almost the 
same. In particular Eq.0} still holds, and the structure 
of the kernel Eq.Q) is not altered, except for a renor- 
malization of £o and Al as shown by Nami&. These two 
length scales are renormalizcd as 



A L -> \ L /VZ (9) 
Co -> £o ■ Z 

where the renormalization factors Z for Pb, Nb and Ta 
arc Z Ph ~ 2.55, Z Nb ~ 2.1, and Z Ta ~ 1.69it 



III. EXPERIMENTAL DETAILS 

The samples investigated in this work were two Pb 
films, a Nb film, and a Ta film. The characteristic prop- 
erties, relevant for the analysis are listed in Table [I] 



Table I: Characteristics of the investigated samples. RRR = 
[R(RT) - R(T C + 1.0)]/R(T C + 1.0) 



sample 


Bcxt (G) 


thickness (nm) 


To (K) 


RRR 






oxide layer 


film 






Pb-I 


88.2 


16(2) 


430(20) 


7.1(1) 


16 


Pb-II 


89.6 


5.8(3) 


1055(50) 


7.21(1) 


23 


Nb 


88.2 


4.2(3) 


310(15) 


9.24(6) 


133 


Ta 


193.4 


1.8(4) 


350(15) 


4.42(2) 


45 



The samples were sputtered directly onto a sapphire 
crystal using 99.999 % pure material. This permits both 
an excellent thermal contact to the cold finger of the 
cryostat and also the application of a bias to the sample, 
which is needed in the LE-/1SR experiments to tune the 
implantation energy of the muons. The Pb films were 
sputtered at room temperature, the base pressure being 
3 • 10~ 7 mbar. The Nb and Ta films were sputtered on 
the substrate maintained at 1000 K. The base pressure 
in the deposition chamber ranged between 2 to 5 • 10~ 8 
mbar. X-ray diffractometry revealed epitaxial growth of 
Nb with a (2 0) growth direction perpendicular to the 
substrate, and the Ta film shows a highly oriented al- 
pha structure most probably with a (2 0) growth di- 
rection. The thickness of the films was determined by a 
high sensitivity surface profiler and Rutherford backscat- 
tering. The surface roughness of these films is given by a 
arithmetic mean roughness value of < 0.7%^. The criti- 
cal temperature T c was measured by means of resistivity 
and susceptibility measurements. The mean free path i 
was estimated according to our resistivity data. 

The oxide layers of Nb and Ta are extremely stable 
as discussed extensively in Refill The diclctric Nb20s 
forms at the surface of Nb and acts as a very good pro- 
tection layer. The typical Nb2 05 layer thickness found 
under the described growing conditions is ~ 5 nm. Ta 
forms a T^Os oxide layer with a typical saturation thick- 
ness of ~ 2 nm, which again acts as a very well protection 
layer for the Ta film. The thickness of the oxide layer as 
determined by our measurements is in excellent agree- 
ment with the findings in the literature. 
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The LE-//SR method makes use of the muon spin rota- 
tion technique 20 - (/iSR) where ~ 100% polarized positive 
muons implanted in a solid sample rapidly thermal- 
ize (~ 10 ps) without noticeable polarization loss. The 
spin evolution of the ensemble after the implantation is 
then measured as a function of time. The evolution can 
be monitored by using the fact that the parity violat- 
ing muon decay is highly anisotropic with the easily de- 
tectable positron emitted preferentially in the direction 
of the fi + spin at the moment of the decay. 

Counting the variation of the decay positron intensity 
N(t) with one or more detectors as a function of time 
after the muon has stopped in the sample, it is possible to 
determine P(t), the time dependence of the polarization 
along the initial muon spin direction. 
The experimentally obtained time histograms have the 
form 




12 3 4 



N{t) = No 



p(ty 

P(0). 



(10) 



iVo is a normalization constant reflecting the total num- 
ber of muons recorded. The exponential describes the de- 
cay of the fj, + and A is the maximum observable asym- 
metry of the decay (theoretically 1/3 in case of 100% 
polarization and when integrating over all positron ener- 
gies). The relevant information about the system under 
consideration is contained in the term Aq P(t)/P(0). 

Unlike conventional /iSR techniques which make use 
of the energetic muons (~ 4 MeV) originating from tt + 
decay at rest ("surface" muons), LE-/1SR makes use of 
cpithermal muons (~ 15 eV) extracted after moderation 
of surface muons from a thin film of a weakly bound van 
der Waals cryosolid (wide band gap insulator By 
re-accelerating the epithermal muons up to 20 keV and 
biasing the sample, it is possible to tune the implantation 
energy in the range of 0.5 to 30 keV and thus to implant 
the muons beneath the surface of any material in a range 
of up to about 300 nm and a spatial resolution down to 
1 nm. Details concerning the fi + stopping distribution 
will be given in Sec lIVI 

The measurements were carried out as follows: A mag- 
netic induction B ex t parallel to the sample surface was 
applied after zero field cooling the sample. The incident 
muon spin was parallel to the surface and perpendicular 

tO £?oxt- 



IV. RESULTS AND ANALYSIS 

The applied magnetic field was determined from the 
Larmor precession frequency of the muon spin at T > T c . 
In this case the measured polarization simply exhibits 
the undamped precession^ of the muon spin ensemble 
P(t) oc cos(7 M B e xti + 0)- In the Meissner state (T < T c ), 
the fx + ensemble, stopping in the surface layer, is sam- 
pling B(z), thus the muons will precess in various fields, 
depending on where they actually stop (distance to the 




f(pS) 

Figure 4: Typical time spectra A P(t, E)/P(0, E). The 
shown data are from sample Pb-II, zero field cooled, T = 
3.05 K, B cxt = 8.82(6) mT. Top: fi + implantation energy 
E = 5.2 keV. Bottom: E = 14.8 keV. 



interface). This yields a damped P(t), reflecting the spa- 
tial distribution of the magnetic field close to the sur- 
face. FigQ] shows some typical time spectra according to 
Eq. HI Up . From these spectra, the magnetic field distribu- 
tion p(B, E) is obtained by Fourier transform 



P(B,E) 



/2vr 



A 



P(0,E) 



cos^Bt + cj>) At (11) 



where 7^ is the gyromagnetic ratio of the /i + and <f> an an- 
gle taking into account the concrete detector geometry. 
E is the implantation energy of the muons. To obtain 
p(B, E) we used a maximum entropy algorithm, which 
proved to be much more robust than the usual Fourier 
transform methods, especially in the case of limited avail- 
able statistica 23 ! 24 ! 25 ! 26 . 

To determine B[z) knowledge of the muon stopping 
distribution n(z,E) is needed. The situation is similar 
to the case of magnetic resonance imaging where a mea- 
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sured nuclear resonance frequency has to be translated 
into a space coordinate. We use the Monte Carlo code 
TRIM. SP2i whose reliability to predict implantation pro- 
files of low energy fj, + in various materials has been con- 
firmed by experimental testing2£. The functional relation 
between n(z, E) and p(B, E) is 



could, in the worst case, only lead to an underestimation 
of present non-local effects or even wash them out. 

The mean values (z) were determined directly from the 
Monte Carlo stopping distribution. The asymmetries, 
measured at various implantation energies, were fitted 
within a Gaussian relaxation model to 



i(z,E) dz = p(B,E) dB 



(12) 



which states nothing else than the probability that a 
given implanted muon with incident energy E stopping 
in the interval [z,z + dz] , will experience a field in the in- 
terval [B, B + dB] with the probability p(B, E) [assuming 
a monotonous B(z)). Integrating Ea. (|12|) on both sides 
yields 



n(C,£)dC 



BO) 



p(P,E)dp 



(13) 



which is, for a chosen z, an equation for B. Since n(z, E) 
can be calculated and p(B, E) can be measured by means 
of LE-/xSR, the magnetic field penetration profile can be 
determined. FigEl shows an example of a B(z) determi- 
nation. 

With this approach it is possible to determine the 
whole B{z) functional dependence from a single mea- 
surement at one specific implantation energy. Still, we 
determined B(z) at various energies which results in a set 
of overlapping curves. This self-consistence check further 
demonstrates the reliability of the Monte Carlo code used 
to determine n(z,E). For a further crosscheck we used 
an additional approach to determine B(z), which allows a 
more rigorous statistical error estimate. With a very nar- 
row stopping distribution, the following two approaches 
also lead to a B(z): (i) Plotting the spatial coordinate 
z p for which the stopping distribution n(z, E) is maximal 
(peak value), against the field value B p where the field 
distribution p(B,E) has its maximum, a B p (z p ) = B(z) 
results for a set of different implantation energies E. (ii) 
Instead of choosing the peak position, the mean values 
(z) = J zn(z,E)dz and (B) = j Bp(B, E) dB can be 
used which again leads to (B) vs. (z) — B(z) for a set of 
different E. However, it has to be taken into account that 
the stopping distribution n(z, E) is not very narrow and 
therefore a more elaborate calculation is needed. The de- 
tails of such an analysis are given in Appendix^ The 
results can be summarized as follows: Both approaches 
introduce minor systematic errors. The systematic er- 
ror of the peak value approach could mimic deviations 
of an exponential decay reminiscent of non-local effects 
{negative initial curvature of log[B(z)]} and therefore 
this approach was excluded. The mean value determi- 
nation is the appropriate choice, since this method only 
produces systematic errors opposing possible non-local ef- 
fects, i.e. an exponential decaying magnetic field profile 
is slightly deformed so that log[_B(z)] has a positive ini- 
tial curvature. The determination based on mean values 



P(t) 
P(0) 



= ^bg exp 
+^sc exp 



~2 <>BGi) 2 

-\ (°sct) 2 



cos^Bcxti 4 
cos(7 M B S ci - 



(14) 



is a phase describing the relative position of the 
positron detectors. The observable asymmetry Aq is 
modeled by 3 contributions A = A sc + A^q + A BS . For 
our present experimental setup Aq — 0.27. Aqc/Aq gives 
the weight of the muons stopping in the superconductor. 
Aqq/Aq is the portion of muons either stopping in the 
oxide layer or in the sample surrounding. These muons 
experience the external field B cxt ■ ^bs I Aq is the portion 
of backscattered muons. The backscattered muons form 
muonium^S (a hydrogen like electron-muon bound state 
with a gyromagnetic ratio 7m u — 1037^). Due to its large 
7Mu, muonium has a much higher precession frequency 
compared to /x + which results, in the present experiment, 
in an instant depolarization (extremely short dephasing 
time, filtering due to sampling of the time signal). This 
leads to an effective reduction in the observable asym- 
metry as written in Ea. ll4fl . The various weights of the 
asymmetries were fixed according to the Monte Carlo 
simulation of n(z,E). Therefore the only free fitting 
parameters are -Bsc, which corresponds to (B) in this 
Gaussian relaxation model approximation, and the two 
depolarization rates <tbg, and asc- FigEI shows a typical 
result together with two B(z) profiles determined from 
Ea. (|13|) for different implantation energies. We would 
like to stress that the mean value approach to determin- 
ing Biz), due to slight systematic errors, rather underes- 
timates the deviation from an exponential magnetic pen- 
etration profile. Still, deviations from the exponential 
behavior of B(z) are clearly present, thus confirming the 
non-local response. 

Additional evidence about the power to detect small 
variations of the magnetic penetration profile is given by 
the previously measured B(z) data of the high temper- 
ature superconductor YBa 2 Cu307_5 (optimally doped; 
for details see Ref l29l) . In this clear cut type II super- 
conductor [£ ~ 1.5 nm, Al = 146(3) nm] a perfect 
exponential B(z) is expected, at least for temperatures 
> 1 K. Below T < 1 K there might be deviations from the 
exponential behavior due to the d-wave pairing charac- 
ter in this compound which can lead to substantial non- 
linear and non-local effectsi^SiSi. We reanalyzed these 
data according to the above discussion. The results for 
T = 20 K are shown in Fig[71 which clearly demonstrates 
the exponential form of B(z). 
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Figure 5: (Color online) Magnetic penetration profile B(z) in sample Pb-II at T = 3.03 K. The implantation energy was 
E = 5.2 keV. Top graph: fi + stopping profile n(z,E) from the Monte Carlo code TRIM.SP. Bottom right graph: p(B,E) 
maximum entropy analysis of Ao P(t) / P(0) [see Eq. llOll l. Bottom left graph: B(z) field determination according to Eg. 1131 . 
The solid line in B(z) is the best fit for the BCS kernel [Eqs.(0J, J7}], whereas the dashed line shows the London limit for 
parameters obtained from the BCS fit. The shaded area of n(z, E) shows the Pb oxide layer. 




20 40 60 80 100 120 140 

z (nm) 



Figure 6: Mean value determination of B(z) (solid dots) in 
sample Pb-II at T = 3.03 K, together with integral approach 
[Eq.O] for E = 5-2 keV (open circles) and E = 14.8 keV 
(open triangles). 
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Figure 7: Magnetic penetration profile B(z) of the high 
temperature superconductor YBa2Cu307_i at T = 20 K 
(T c = 87.5 K). The B(z) obtained from the mean values is 
displayed by solid circles whereas the other points show the 
integral determined field profile based on Ea. l)13^ . 



After having established the reliability of our approach 
to determine B(z), we turn to the discussion of the data. 
FigEl shows a collection of our results for Pb; FigllUI 
shows the Nb and Fig^] the Ta data. The low tem- 



perature data of both Pb samples show a clear deviation 
from an exponential decay law, indubitable evidence for 
the presence of non-local effects. Furthermore, increas- 
ing the temperature leads to less-pronounced curvature 
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as expected, since very close to T c non-local effects should 
disappear altogether (strong temperature dependence of 
A for T — > T c ; weak one of £ [see FigQ]). Unfortunately, 
the range where a sign reversal of B(z), predicted by the- 
ory, should appear, is experimentally not accessible yet. 
To compare with theory, the data were fitted according 
to Eq.0 for the Pippard kernel [Eq.@] and the BCS 
kernel [Eqs.Q,©], respectively. Since the temperature 
dependence of A is undefined within the Pippard model , 
we chose the two- fluid approximation A(i) = Ao/vl — t 4 , 
with t — T/T c . X is an effective London penetration 
depth, taking into account corrections due to the scatter- 
ing of electrons. In the clean limit (£ — > oo) Ao = Al- The 
temperature dependence of A within the BCS model is 
given implicitly by the Eqs.|Q),ljBJ), where again Ao is used 
to indicate an effective London penetration depth. The 
fits to the data are shown in Fig[S]and were obtained in 
the following way: for the low-temperature high-statistic 
data it was possible to fit Ao and £o- For the rest of the 
data, £o w &s fixed to the value found at low tempera- 
ture. This was necessary, since by approaching T c , the 
magnetic penetration profile becomes more and more ex- 
ponential due to the strong temperature dependence of A. 
As a consequence, the ^-dependence of B(z) weakens re- 
sulting in a drastically growing uncertainty to determine 
the parameter. 

In the models used to analyze B(z) a temperature de- 
pendence of the magnetic penetration A(T) is assumed 
(Pippard Sec III A|) or implicitly given (BCS Sec lIIBjl . 
If these models describe the temperature dependence of 
A(T) correctly, Ao should be a constant value within the 
error bars for all the data sets. FigOshows a graph where 
Ao is plotted versus the reduced temperature t = T/T c for 
both models. For the two-fluid temperature dependence, 
assumed in our Pippard model, Ao is indeed tempera- 
ture independent, except for data close to T c . The weak 
coupling BCS temperature dependence for the magnetic 
penetration length, however, produces a clear temper- 
ature dependent Ao meaning that the real temperature 
dependence of the energy gap A(T) does not follow ex- 
actly the weak coupling prediction. Tabic [H] shows the 
collected results. 

In the case of Nb (Fig llOfl the limited statistics of the 
data do not allow us to fit £o directly. For the analysis it 
was therefore fixed to its literature value (see Table [n}. 
This system is at the borderline between the non-local 
and local regimes. Therefore B(z) was also analyzed with 
a simple exponential model, leading to Ap Xp = 33(3) (nm) 
which is similar to the non-local value given in Table ITTI 
The x 2 suggests that also for Nb the non-local regime 
is the appropriate one, but definitely higher quality data 
are needed to resolve this issue. 

The Ta data (Fig lll(l show a more pronounced devi- 
ation from an exponential decay law. As for Nb it was 
necessary in this case to fix £q to the literature value. 
Table ITTI summarizes the results. In the data analysis we 
also considered the possibility of spurious effects related 
to surface roughness or thickness variation of the film 



Table II: Results of the analysis for Pb, Nb and Ta. The 
values for Aq CS,p given, are the low temperature values. The 
"dead layer" corresponds to a shift of the ^-origin of B(z). 



sample 


compound 


\BCS,P 
A 

(nm) 


Co 
(nm) 


dead layer 
(nm) 


Pb-I 


Pb 


59(3) 


90(5) 


6(1) 


Pb-II 


Pb 


55(1) 


90(5) 


3(2) 


Nb 


Nb 


27(3) 


39(fixed)22 


2(2) 


Ta 


Ta 


52(2) 


92(fixed)2£ 


3(1) 



due to the presence of terrace-like structures. If these 
structures have a lateral size smaller or comparable to 
the coherence length, the screening current will be very 
ineffective, thus resulting in a dead layer at the surface. 
We effectively find a better fit to the data assuming the 
presence of such a dead layer (see FigH^l and Table ITTft 
in addition to the previously mentioned oxide layer. The 
other case of surface terraces with lateral extension ^> £o 
would mimic a thickness distribution of the film. Numer- 
ical simulations have shown that such a thickness varia- 
tion would only marginally affect the magnetic penetra- 
tion depth. A closer look at the Ta curves shows that 
compared to the other systems the fit based on the BCS 
or Pippard model does reproduce the data in a less sat- 
isfactory way. One possible origin for this discrepancy 
is the neglect in both models of the suppression of the 
supercurrent density on approaching the surface. Fig 1 121 
sketches the situation. The models discussed in SeclTllas- 
sume ns to be constant up to the surface (dashed line). 
However, from the Ginzburg-Landau theory it is known 
(see e.g. Ref l32l) that close to the surface the order pa- 
rameter, and hence ns, are reduced. Such a reduction 
in ns makes the screening less effective, which, on top of 
the non-local effects, would lead to a more pronounced 
deviation from the exponential penetration profile than 
estimated by Eq.QJ or its diffuse scattering counterpart. 
In the presently investigated samples we expect such an 
effect to be more pronounced in the cleaner supercon- 
ducting Ta than in Pb (where also small deviations from 
the theoretical curves are possibly present in Fig|SJ be- 
cause of its lower effective k. In order to quantify these 
effects in a comparison with our experimental findings, it 
would be very useful to have a self-consistent theoretical 
description of the physics outlined above. 



V. SUMMARY 

We have performed low energy muon spin rotation 
spectroscopy experiments (LE-/iSR) on Pb, Nb and Ta 
films in the Meissner state. The magnetic penetration 
profile B(z) into the superconductor has been deter- 
mined on the nano-meter scale, thus providing a model 
independent measure of the Meissner screening profile. 
B(z) shows clear deviations from the simple exponen- 
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Figure 8: (Color online) Magnetic penetration profiles for Pb at various temperatures. Left graph: Data from sample Pb-I. 
Right graph: Data from sample Pb-II. The solid lines are BCS fits to the data, whereas the dashed line represents B(z) = 
Bext exp( — z/X), where the X from the BCS fit is used. 



o 




Figure 9: Ao determined at various temperatures t = T/T c 
for the Pb data. The temperature dependence of A in the 
Pippard model was chosen as 1/A 2 oc \/l — t 4 , whereas in the 
BCS model the temperature dependence is given implicitly 
by Eqs.@,(|7J,(|SJl. The curves are guides to the eyes. 



tial decay law as expected in the Pippard regime where 
the relation between the current density and the vec- 
tor potential is non-local in nature. Analyzing the data 
within the Pippard and BCS model, the coherence length 
£o and the London penetration length Al could be de- 
duced. Furthermore we could show that the two-fluid 



model approximation for the temperature dependence of 
A oc 1/ y/l — (T/T c ) 4 is closer to the experiment than 
the BCS one which is A oc l/y/K B cs(q ->0,T,£^ oo). 
These experiments based on the local profiling of the 
magnetic field penetrating at the surface are the first 
measurements showing the non-local nature of super- 
conductivity in Pb, Nb and Ta on the nanometer scale, 
and verify directly and quantitatively the longstanding 
predictions of Pippard. B(z) for Ta deviates substan- 
tially from the simplest theoretical approach, which we 
attribute to the neglect of the reduced superfluid density 
on approaching the surface. 
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Figure 12: Effect on B(z) of the suppression of ns(z) at the 
surface. The shaded area represent the sample cross-section. 
The left side shows a surface with a terraces of typical size 
< £o , which results in a "dead layer" . The dashed lines show 
ideal superfluid density ns(z) and magnetic penetration pro- 
file B(z), respectively. The solid lines are the corresponding 
more realistic ones. 



Figure 10: Mean value determination (solid dots) in Nb at 
T = 2.95 K and T = 7.6 K, together with the integral ap- 
proach. The solid lines show BCS fits. 




Figure 11: Magnetic penetration profiles for Ta at various 
temperatures. The solid lines show BCS fits. 
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Figure 13: Model stopping distribution n(z, E) of low energy 
muons. 



Appendix A: MODELING THE 
DETERMINATION OF THE MAGNETIC FIELD 
PROFILE 

As pointed out in the paper the integral determination 
of B(z) [Eq. l|13|l ] is the most effective method since it is 
able to generate a complete curve from a single implan- 
tation energy. As a cross-check of the use of the maxi- 
mum entropy method to determine p(B, E) and B(z) via 
Eci. (|13[) . where the determination of the statistical error 
is an unresolved issued, we analyzed the /iSR data di- 
rectly in the time domain which permits estimation of 
the statistical errors. The time domain modeling used in 
our analysis has the disadvantage of introducing system- 
atic errors and only one of the discussed approaches is 
acceptable for the investigation of non-local effects. 

Two sets of models, which will be discussed separately 
in the following paragraphs, can be treated analytically. 
The resulting trend is similar in both cases, meaning 



that the "mean value determination" fSec IA Q slightly 
overestimates the real B(z), whereas the "peak value 
determination" f Sec I A 2\i always slightly underestimates 
B{z). In the situation where one is looking for non- 
local effects, characterized by an initial negative curva- 
ture of log \B{z) I 'Bextl, the mean value reconstruction is 
the method of choice because only it can produce system- 
atic errors with a positive curvature, and henceforth, in 
the worst case, can only hide the presence of a negative 
curvature. All these results are compiled in Figll4l 



1. Mean Value Determination of B(z) 

The muon stopping distribution can only be calculated 
by Monte Carlo codes. For the following discussion a 
simple stopping distribution mimicking a realistic one is 
needed. We take as a model stopping distribution n(z, E) 

n(z,E) — n Q (z - z) z 4 , ze[0,z o ], (Al) 
where no = 30 Zq 6 is the normalization factor and zq 
the maximum distance, which an implanted muon can 
reach. This function is an acceptable approximation for 
a realistic implantation profile fFigl!3[l. 
The mean value of z is 

(*> = f*Q. (A2) 

Assuming further 

B{z) = B cxt exp(-z/A) (A3) 
one gets, utilizing the identity n(z, E) dz — p(B, E) dB, 



(B) 



dzB(z)n(z,E) 

= B CKt n X 2 [24A 3 (z - 5A) + e~ Zajx {z\ + 8Xz^ + 36A 2 z^ + 96A 3 z + 120A 4 } 
= £ cxt cxp(-(z) /A) + /((*)), /««))> 

I 



with lim f((z)) = and f((z)) < (B)((z)) in the first 

(z>->0 

two decades as shown in Fig^J Notice that the curva- 
ture of f((z)) is positive. 



In order to have an estimate not only for an exponential 
B(z), also the following B(z) (resembling a non-local field 
profile close to the solution in the extreme anomalous 



limit 1 - ^ ^> A) was analyzed 
B(z) = B cxt exp 



V3, 
2 



Qz 



where Q 



3tt 1 



4 A^o 



1/3 



(A4) 



(A5) 



Also here an exact solution 



for the mean value can be given 
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where the coefficients Ci are 



(B)/B ext = co + ciexp \-~l[^-Q( z )J C0S [j^Q{ z ) ) + 

+c 2 cxp Q(z) \ sin (L Q( z ) ) (A6) 

J 



10 ^ 1 ' / v io 



co = - 



Cl 



C2 



1125000 (7V3Q(z) -50) 

117649 (Q(z)) 6 
375 [150000 + 84000 V3Q(z) + 44100 (Q(z)) 2 



2401 (Q(z)) 4 



117649 {Q(z)Y 

375 [12000 + 6300%/3 Q(z) + 3920 (Q(z)) 2 + 343 V / 3(Q(^)) 3 ] 



16807 (Q<») 5 
I 



These formulae look rather ugly but as can be seen in 
Fig ll4l the result is very similar to the simpler case dis- 
cussed before. 



Peak Value Determination of B(z) 



with Bo the position of the Gaussian peak. SB is the 
width of the distribution. Assuming furthermore an ex- 
ponential B(z) [Ea. (|A3|l ]. one arrives at the following 
stopping distribution 



First lets assume the model stopping distribution from 
Ea. (jAl|l . The peak value of the spatial coordinate is 



z p = - z Q . (A7) 
The corresponding peak position is obtained from 



n(z,E) 



1 B Q 



2n\ SB 



exp 



1 { B cxt e- Z / X - B \ _£ 
SB J A 
(A10) 



The peak values are than given by 



p(B,E) = n(z,E) 
A 



dB 



dz 



B 



(A \n{B ext /B)f (z - A ln(B ext /B)) . 



with B € [B cy± e~ z °/ X , B cxt \. Since we are only interested 
in the peak position, the normalization factor was sup- 
pressed. It follows that 



Bp = S oxt exp 



gT- ( z o - 5A + V(^o + 5A) 2 - 4z 



B cxt exp(-z p /A) - g(z p ), 



(A8) 



with g(zp) > and lim Zp ^o g{ z p) = 0. The results found 
so far are compiled in Figll4l 

Another ansatz, which is extremely close to the our 
analysis is the following: Instead of assuming a model 
distribution n(z, E) [Ea. (|Al|) ]. one assumes that p(B, E) 
is Gaussian, i.e. the time dependent LE-/1SR signal is 
Gaussian damped. 



p(B,E) 



1 



2irSB 



exp 



1 ( B-Bp 

2 V SB 



(A9) 



Bp — B a 



A In 



B c 



2SB 2 



which leads to 



Bp(zp) — B cxt e" 



z p /\ 



1 - 



Bo + iSB - Bo 



SB 
B rx t 



+2z p /\ 



, (All) 



showing the same trend as in the previously discussed 
case. 

In conclusion, one finds that the peak value approach 
to determine B(z) leads to systematic errors that can 
mimic a magnetic penetration profile of a non-local su- 
perconductor and therefore has to be excluded from such 
an analysis. The mean value determination of B(z) also 
introduces systematic errors, however they tend to dimin- 
ish real non-local effects and are therefore less dangerous. 
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Figure 14: (Color online) Left: Analytic results for mean- and peak value determination assuming an exponential B(z). 
Right: Extreme anomalous limit result for the mean value determination of B(z). 
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